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Problem Statement 


Spatially map sound with 
qualitative/quantitative accuracy 

Microphone arrays allow for spatial 
separation of distinct sources 

Many existing processing methods 




Downstream view of Hybrid 
Wing Body model inverted on 
test stand with phased 
microphone array overhead 
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Existing Spatial Sound Mapping Methods in Aeroacoustics 


• Cross-spectral 

o Beamforming 
o CLEAN-PSF [Hogbom 1974] 

o Spectral Estimation Method [Blocodon & Elios 2003] 
o DAM AS and DAMAS-C [Brooks & Humphreys 2004, 2006] 
o CLEAN -SC [Sijtsmo 2007] 

o Functional Beamforming [Dougherty 2014] 

• Eigenspace 

o Generalized Inverse [Suzuki 2008] 
o Orthogonal Beamforming [Sarradj 2010] 

• Wavespace 

o Wavespace Deconvolution [Bohr& Cottofesto 2012] 



Research Objectives 


Use: 

>Single-channel filtering based on user-defined spatial constraints 


As: 


Building-block for existing spatial mapping techniques 


■ For 

>Qualitative/quantitative improvement in accuracy 



Microphone Array Filtering 

• Goal: Filter channel data for more accurate spatial sound estimation 


Filter dashed-red to 
better estimate solid-white 





Modified Wiener-Hopf Eg. for Spatial Filtering 

• Combine : Optimal, least-squares filtering of Wiener-Hopf Eq. + Spatial filtering 

• Only constraint = User-defined phase 

• Modified weight vector -> Filters + Prevents targeted signal cancellation 
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*Constrained Spectral Conditioning (CSC)* 


• Conditioned Spectral Analysis 1 with modified Wiener-Hopf Eq. (WB) becomes CSC 2 

• Optimal, spatially-constrained, least-squares filtering for Fourier Transforms of microphone outputs 


Wiener-Hopf Eq. 


Modified Wiener-Hopf 
Eq. via Orthogonal 
Projection 


Constrained 

Spectral 

Conditioning 
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1 Bendat & Piersol 1986, 2 Spalt 2014 
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Basic Spatial Sound Mapping 



Frequency-Domain 
Beamforming (FDBF) 

(a.k.a. Delay-and-Sum (DAS)) 



FDBF via CSC 


CSM formed 
from CSC 
outputs = 
"CSC-CSM" 


FDBF via 
CSC 
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FDBF B eamw i dth and Sidelobes: Simulate d Point Source 

• Point source "measured" with SADA 1 : 60" from array face,/ = 10 kHz, SNR = 20 dB 




Humphreys et al. 1998 
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FDBF Beamwidth & Highest Sidelobe Level vs. Frequency 

• CSC performance dependent on SNR & frequency 




SNR, 

dB 

Highest Sidelobe of FDBF via CSC- 
CSM relative to highest sidelobe of 
FDBF via initial CSM, dB 
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CSC Observations 


• Constrained Spectral Conditioning (CSC) 

• Single-channel processing -> Building-block for existing, algorithms as it 
processes the Fourier Transforms of microphone outputs 

• Uses only relative phase differences as constraints 

• CSC success dependent on: 

• Frequency 

• Source field 

• Solid angle 

• Microphone layout 

• Signal-to-Noise Ratio (SNR) 



CSC Observations, cont. 


• CSC output datasets ("CSC-CSMs") are estimates 


• FDBF using CSC-CSMs 

• Non-integrateable (due to inaccuracies) 

• Non-linear Cannot be deconvolved "easily" 


r Modified approach needed for accurate spatial sound mapping 



Iterate 


Advanced Spatial Sound Mapping via CSM Decomposition 

• CLEAN-PSF 1 





PSF 
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CLEAN-PSF 




max 


i-1 


FDBF to locate 
max location 
on map 


Subtract CSM 
estimate 


Iteration stop 
criterion 


Output 

source 

map 


Pros 

o Takes advantage of 
"uncovered" information 


Cons 

o No sidelobe discrimination 
o Relies on the PSF for source 
propagation model 


1 Hogbom 1974 
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Iterate 


Advanced Spatial Sound Mapping via CSM Decomposition, cont. 

• CLEAN-CSC 1 



’'Hip)'-' 8 '** 


FDBF to locate 

ses 

max location 
on map 


CSC-CSM 


G i+1 = &- <pYg csc 

s ma x£SC 

Subtract CSM 
estimate 

Stop if : ^ G l > 

Z 6 '" 1 

Iteration stop 
criterion 

i 

Y CLEAN -CSC = <P / J 

ft esc 

Output 

source 

i-1 


map 


Pros 

o Higher location/level accuracy than 
FDBF 

o Sidelobe discrimination 
o Does not use PSF magnitudes 

Cons 

o CSC-CSMs have inaccuracies 
o PSF phase still used 
o Cannot take advantage of 
"uncovered" information 


1 Spa It 2014 
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Iterate 


Advanced Spatial Sound Mapping via CSM Decomposition, cont. 


• CLEAN-SC 1 
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CLEAN-SC-CSM 
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Y CLEAN SC ~ <P ^ 
i=l 


FDBF to locate 
max location 
on map 

Form CSM 
estimate 


Subtract CSM 
estimate 

Iteration stop 
criterion 

Output 

source 

map 


Pros 

o Adaptively defines CSM estimate 
magnitudes/phases 
o Takes advantage of "uncovered" 
information 

Cons 

o No sidelobe discrimination 
o Multiple/distributed sources bias 
CSM estimate 

o Stronger sources bias weaker 
sources 

o Inaccurate for coherent sources 
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1 Sijtsma 2007 



CLEAN-CSC + CLEAN-SC ± CLEAN-CSC-SC 1 

• CSC -> FDBF and CSM at max locations 

• CLEAN-SC -> Further decompose CSC-CSM -> CSC-SC-CSM 

• Improves CSC-CSM magnitude estimates 

• Corrects deviations in initial phase definitions 

• CLEAN-SC only used once original CSM is sufficiently decomposed 

• Improves dynamic range 



r 
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CLEAN-SC 


SC-CSM 


Initial CSM — I CSC-CSM | < ICSMI ? 


H 


Yes 

L CLEAN-CSC 


CSC-CSM - CLEAN-SC 


CSC-SC-CSM 


*S palt et al. 2015 
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Simulation: Modified PSF, Incoherent Sources 


467 point sources "measured" with JEDA 1 : 72" from array face,/= 15 kHz, SNR = 20 dB 

True Source Distribution CLEAN-PSF 
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Simulation: Modified PSF, Incoherent Sources, cont. 


CSM usage during decomposition 



CLEAN-SC vs. CLEAN-CSC-SC 
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Preliminary Jet Noise Results 



JEDA (left) positioned at 90° with respect to the 
jet exit plane in the JNL. 


• Single-stream, convergent nozzle 

• Exit diameter 2.67" 

• Supersonic, cold jet at Mach 1.48 

• Wind tunnel co-flow at Mach 0.1 

• Array ("JEDA") 6 ft from jet 
centerline 
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CLEAN-CSC-SC CLEAN-SC 


Preliminary Jet Noise Results, cont. 

* No CSM Diagonal Removal, No CSM Weighting 



8 kHz 1/3 oct. Mean Array Single Mic level = 113.1 dB 
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CSC Conclusions 


• Building block for existing algorithms 

• Improves result accuracy under incoherent source conditions 

• Not "plug-and-play" for advanced spatial mapping algorithms 



CLEAN-CSC-SC Conclusions 


• More qualitatively/quantitatively accurate than CLEAN-SC for 
incoherent sources 

• More analysis needed when source coherence exists 
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CLEAN-CSC-SC 1 


• Spatial sound 
mapping via 
decomposition of 
the initial CSM 
using CSC 2 and 
CLEAN-SC 3 




s max ,CSC 


h- 


^ 5 JJU IX 


C i max £scw. 


1 1 + V »4 HWt 

\J *mnx s tnaxy 


S,nnx t-J' "Imai 

y t max ,csc 


G Smax ,csc-sc - ^ n ^sc[(fc^csc-5c)(Af /:w-sc)] 


if \Gt mBX xsc\ < |£*| 


G l+l = G‘ - (pG; csc-SC 


else 


g j+1 = &- ^n_[(%,^c)(^_5c)] 


1 Spa It et al. 2015 
2 Spalt 2014 
3 Sijtsma 2007 
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1. Beamform to locate max grid location 

2. Calculate CSC beamform estimate 

at max location using a CSC-CSM 

3. Form normalized steering vectors from CSC 

beamform and CSC-CSM using CLEAN-SC 

4. Calculate CSC-SC-CSM at max location using 

normalized steering vectors 

5. Lower CSM energy is deemed more 

accurate 


6. Update decomposed CSM accordingly 


7. Stop if CSM energy increases or remains 
unchanged 
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Outline 


1. Introduction 

2. Research Methodology 

3. Simulated Data Analysis 

4. Experimental Data Results 

5. Contributions 


6. Future Work 
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Research Methodology 

CSC Extension to Full Array 

• *Optimum reference channel for use in arrays*: 

1. Maximize undesired signal cancellation 

2. Prevent amplification of noise 



m' 0 (m,s ) = max\WB\ 
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m! ^ m 
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Research Methodology 

*CSC Iterative Processing Algorithm * 
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• *Stop if*: 

1. Processed channers magnitude > channel's initial magnitude 

2. Coherent signal between channels < noise floor between channels 
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Single-channel array inputs 


Undesired Signal Cancellation via Spatial Filtering 



• Generalized Sidelobe Canceller (GSC ) 1 ' 2 

• Filtering method designed to attenuate all signal except that from a user-defined point in space/direction 

• Filtering performed on the synthesized array data 



targeted spatial 
location 


Adaptive filter which modifies 
x B in order to maximize y B 


Output beamform 


1 Hanson & Lawson 1969, Griffiths & Jim 1982 
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Single-channel array inputs 


Undesired Signal Cancellation via Spatial Filtering, cont. 


• Modified GSC 1 



Output 

beamform 


Blocking/filtering performed at the s/ng/e-cfranne/ level 


1 Liu & Van Veen 1991 
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Background Subtraction" for CSC* 
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